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Abstract 
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I. INTRODUCTION 



It is well know that warm inflation, as opposed to the conventional cool inflation, presents 

n 

the attractive feature that it avoids the reheating period ]1J. In these kind of models dissipa- 
tive effects are important during the inflationary period, so that radiation production occurs 
concurrently together with the inflationary expansion. If the radiation field is in a highly 
excited state during inflation, and this has a strong damping effect on the inflaton dynamics, 
then, it is found a strong regimen of warm inflation. Also, the dissipating effect arises from a 
friction term which describes the processes of the scalar field dissipating into a thermal bath 
via its interaction with other fields. Warm inflation shows how thermal fluctuations during 
inflation may play a dominant role in producing the initial fluctuations necessary for Large- 
Scale Structure (LSS) formation. In these kind of models the density fluctuations arise from 
thermal rather than quantum fluctuations [2]. These fluctuations have their origin in the 
hot radiation and influence the inflaton through a friction term in the equation of motion 
of the inflaton scalar field jsj]. Among the most attractive features of these models, warm 
inflation end when the universe heats up to become radiation domination; at this epoch the 

n 

universe stops inflating and "smoothly" enters in a radiation dominated Big-Bang phase pj. 
The matter components of the universe are created by the decay of either the remaining 
inflationary field or the dominant radiation field 4|. 

In the Dvali-Gabadadze-Porrati (DGP) model 5j the induced gravity brane-world was put 
forward as an alternative to the Randall-Sundrum (RS) one-brane model in which general 
relativity was recovered, also despite an infinite extra dimension, but without warping in 
5-dimensional Minkowski space-time. In the DGP model, the gravitational behaviors on 
the brane are commanded by the competition between the 5D curvature scalar in the bulk 
and the 4D curvature scalar on the brane. In contrast to the RS case with high energy 
modifications to general relativity, the DGP brane produced a low energy modification. In 
the DGP model, according to the embedding of the brane in the bulk, there appear two 
branches of background solutions. For a review of the phenomenology of DGP model, see [3] 
and inflation models and reheating in this scenario were studied in j^l, 0, 10. 11]. 

Usually, in any models of warm inflation, the scalar field, which drives inflation, it is the 
standard inflaton field. As far as we know, a model in which warm inflation on a warped 
DGP brane has not been yet studied. The main goal of the present work is to investigate 
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the possible realization of a warm inflationary universe model, where the energy densities 
(inflaton-radiation) are confined to the brane in DGP model. In this way, we study warm- 
DGP model and the cosmological perturbations, which are expressed in term of different 
parameters appearing in our model. These parameters are constrained from the WMAP 
three year data [121 ] . 

The outline of the paper is a follows. The next section presents a short review in the 
Friedmann equation on the warped DGP inflation model. In Section II I II we present the 
warm inflationary phase on the DGP brane. Section IIVI deals with the scalar and tensor 
perturbations, respectively. In Section |V] we use an exponential potential for obtaining 
explicit expression for our models. Finally, See l VII summarizes our findings. We chose units 
so that c — h—1. 



II. FRIEDMANN EQUATION ON THE WARPED DGP BRANE 



We start by writing down the Friedmann equation on the warped DGP brane, by using 
the Friedmann-Robertson- Walker (FRW) metric. It becomes 

k 1 



p + Po(l + eA(p,a)) 



(1) 



a" 3/i 2 

where H is the Hubble parameter, a represents the scale factor, p is the total energy density, 
and k is the constant curvature of the three-space of the FRW metric. The p parameters 
denotes the strength of the induced gravity term on the brane. Also, the e parameter 
becomes either +1 or —1 representing the two branches of this model. For e = — 1 the brane 
tension can be assume to be positive, while for e = +1 this is negative. A is defined by 



A 



where 

A = 

and A is defined by 



2r] 



p 2 A 
Po 



p = m x + 6- 



,mi 



A=^( (5) A + ^A 2 ), 



6m® 
poP 2 



(0 < 7/ < 1), 



(2) 



(3) 



(4) 



where k$ is the 5-dimensional Newton constant, ^ 5 ^A is the 5-dimensional cosmological con- 
stant in the bulk, and A is the brane tension. Note that there are three mass scales, p, 
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A 1 / 4 and 



m 5 



-2/3 



Also, the quantities Sq is a constant related to Weyl radiation. 



Since we are interested in the inflationary dynamics of the model, we will neglect the cur- 
vature term and the dark radiation term. We shall restrict ourselves to the RS critical case, 
i.e. A = 0. Then Eq.flTJ becomes 

2 VPm/2 



1 



p + p + epo(l + — ) 
Po 



(5) 



Note that in the ultra high energy limit where p ^> po m\, Eq.(|5]) results to be 

1 



3/i 2 



(6) 



In the intermediate energy region where p <C po but p ^> m\, for the branch with e = — 1, 



the Friedamnn equation reads H 2 



L8mf V P 2m\ 



_ V 2 p 2 
6m| " 4m| " 



m 4 <C po i.e. in low energy limit, Eq.fJSJ) becomes 

1 



3p2 



P + - 

Po 



p 



Finally, when p <C 



(7) 



where p p is the effective 4-dimensional Planck mass and is given by p 2 , = p? j (1 — if). 

In the following we will consider a total energy density p = p^ + p 7 where <fi corresponds 
to a self-interacting scalar field with energy density, p^, given by, p^ = |0 2 + V((f>) and p 7 
represents the radiation energy density. 



III. WARM-DGP INFLATIONARY PHASE 



The dynamics of the cosmological model in the warm-DGP inflationary scenario is de- 
scribed by the equations 

4>+ 3Hj> + V i4> = -r 0, (8) 

and 

p 7 + AH Pl = (9) 

Here T is the dissipation coefficient and it is responsible of the decay of the scalar field into 
radiation during the inflationary era. T can be assumed to be a constant or a function of 
the scalar field or the temperature T or both Here, we will take T to be a function 
of 4> only. At the end subsection B of section V, we briefly describe the situation in which 
T = const. In the near future we hope to study more realistic models in which T not only 
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depends on but also on T, expression which could be derived from first principles via 
Quantum Field Theory approach [3, Q- On the other hand, T must satisfies T = /(</>) > 
by the Second Law of Thermodynamics. Dots mean derivatives with respect to time and 
V A = dV^/dcj). 

During the inflationary epoch the energy density associated to the scalar field is of the 
order of the potential, i.e. p^ ~ V, and dominates over the energy density associated to the 
radiation field, i.e. p^ > p 1 . Assuming the set of slow-roll conditions, i.e. <f> 2 -C V(<p), and 
4> <C (3H + T)0 the Friedmann equation ([5]) reduces to 



H 2 



1 



V + p + ep (A 2 + 2 -^f 2 



and Eq. ([8]) becomes 



3H \ 1 + r 



-V, 



where r is the rate defined as 



r 

3H' 



(10) 



(11) 



(12) 



For the high (weak) dissipation regimen, we have r ^> 1 (r < 1). 

We also consider that during warm inflation the radiation production is quasi-stable, i.e. 
p 7 <C AH pry and p 7 <C Tip 2 . From Eq.Q we obtained that the energy density of the radiation 
field becomes 

T0 2 



AH 



(13) 



which could be written as p 7 = <jT^, where a is the Stefan-Boltzmann constant and T r is 
the temperature of the thermal bath. By using Eqs. ffTT]) . ffT2l and ffTB"]) we get 



p y = a 



r p 



4 (1 + r) : 



V": 



V + Po + epMl + ^Y 1 



(14) 



Introducing the dimensionless slow-roll parameter we get 



H 

H 2 



2 1 + r 



V 



1 + er](A 2 + 2r) V/po) 



-1/2 



l + f(l 



e(A 2 + 2 V ^y/ 2 



(15) 



and the second slow-roll parameter a becomes 

H V., 



a 



till 3(1 + r)H 2 ' 



(16) 
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We see that for r = (or T = 0), the parameters e and a given by Eqs. (ll5l) and (TT6]) 
respectively, a re .educed to the typtca, e„n for cool inflation in the DGP brane^ 
Note that the term in the bracket of Eq. (1151) is the correction to the standard warm infla- 
tionary model. Also, this term loses the condition for inflation, i.e. e <C 1, and when e = 1 
and tightens the condition when e = —1, in contrast to the RS model where the correction 
term always loses the inflationary condition. 

It is possible to find a relation between the energy densities p 1 and p^ given by 



Pi 



2(1 + r) 



P<t> + p (l + e[A 2 + 277 



>S1 1/2 ) 



l + e V (A 2 + 2 v /Vpo)- 1/2 



2(1 + r) 



V + p (l + e[Al + 2r ] ^) 



l + er l (A 2 + 2r ] V/p )- 1 / 2 



Recall that during inflation the energy density of the scalar field becomes dominated by the 
potential energy, i.e. p^ ~ V . 

The condition which the warm inflation epoch on a DGP Brane could take place can be 
summarized with the parameter e satisfying the inequality e < 1. This condition is analogue 
to the requirement that a > 0. The condition given above is rewritten in terms of the 
densities by using p 7 , we get 



+ p (l + e[^ + 27 7 ^l 1 ^ 



P0 J 



1 + eq(Al + 2r] p^/po 



-1/2 



2(1 + r) 
> Pi- 



(17) 



Inflation ends when the universe heats up at a time when e ~ 1, which implies 



Ya 

V 



1 2 



1 + erj{Al + 2?7 V/ Po ) 



-1/2 



1 + ^(1 + 6(^ + 2^)1/2) 



;i + r). 



(18) 



The number of e-folds at the end of inflation is given by 

-<Pf h 2 



N 



r<pf h 2 
-3/ — (l+r)#'. 



(19) 



In the following, the subscripts * and / are used to denote to the epoch when the cos- 
mological scales exit the horizon and the end of inflation, respectively. 

IV. PERTURBATIONS 



In this section we will study the scalar and tensor perturbations for our model. Note 
that in the case of scalar perturbations the scalar and the radiation fields are interacting. 
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Therefore, isocurvature (or entropy) perturbations are generated besides of the adiabatic 
ones. This occurs because warm inflation can be considered as an inflationary model with 
two basics fields Q]- In this context dissipative effects can produce a variety of spectral, 



ranging between red and blue 



(2, UM, and thus producing the running blue to red spectral 

n 



suggested by WMAP three-year data[12l|. 

As argued in Ref.[10|] for the DGP brane (see also 



could be written as Sh = I 1 

O (p 



f 50 0. 

From Eqs. tfTT!) and (|T2|) . the latter equation becomes 

36 H*r 2 „ , 



the density perturbation 



8 2 



H 



25 V 2 



(20) 



The scalar field presents fluctuations which are due to the interaction between the scalar 
and the radiation fields. In the case of high dissipation, the dissipation coefficient T is much 
greater that the rate expansion H , i.e. r = T/3H ^> 1 and following Taylor and Berera[l^], 
we can write 

where the wave-number kp is defined by kp = \/TH /V = H y/3r > H, and corresponds to 
the freeze-out scale at which dissipation damps out to the thermally excited fluctuations. 
The freeze-out wave-number kp is defined at the point where the inequality V^^ < TH, is 
satisfied 



From Eqs. ([20} and (J2H) it follows that 

r9 18^3 



'H 



25 7T 2 



H 5 T r r 5 ' 2 



V 2 



(22) 



The scalar spectral index n s is given by n s — 1 



d In 8^ 

, H , where the interval in wave number 

dink ' 



is related to the number of e-folds by the relation dlnk((j)) 
get 



-dN((p). From Eq.([22]), we 



1 - 



he - 2a + 



hr A V A 
6r 2 H 2 



i-5£ 

2 



25 



,4> L ,4> 



2HT 2 



(23) 



where, the slow- roll parameters e and a, (for r> 1) are given by 



EL 

2r 



V 



l + er ] (A 2 + 2r ] V/p c 



-1/2 



^(l + 6(^ + 2^)V2) 



and a 



V, 



3rH 2 



(24) 



respectively. 
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One of the interesting features of the three-year data set from WMAP is that it hints 
at a significant running in the scalar spectral index dn s /d\nk = a s 12]. Dissipative effects 
themselves can produce a rich variety of spectra ranging between red and blue [2j, [l5j. From 
Eq. ([2"3l we obtain that the running of the scalar spectral index becomes 



25 



(lnr) 



3V, 



2 P-^-2w{^- 2 f)\- (25) 
In models with only scalar fluctuations the marginalized value for the derivative of the 
spectral index is approximately —0.05 from WMAP-three year data only [l^ . 

As it was mentioned in Ref. [201] the generation of tensor perturbations during inflation 
would produce stimulated emission in the thermal background of gravitational wave. This 
process changes the power spectrum of the tensor modes by an extra temperature depen- 
dently factor given by coth(/c/2T). The corresponding spectrum becomes 

2 



16tt fH\ 
At = — — coth 
9 fi 2 \2n J 



' k ' 


4 


_2T_ 





V + p + ep (A 2 + 



where the spectral index n g , results to be given by n 



d\nk 



Po 



In 



/2 





" k ' 


coth 


2T_ 



coth[fc/2T] 



(26) 



—2 e. Here, 



we have used that A 2 g oc k n <> coth[fc/2T] [20j. 

For r> 1 and from expressions ( 1221) and ( 1261) we may write the tensor-scalar ratio as 

~_g?LJ V A coth^l 



(27) 



1 /2 

Here, 5h = 2 /5 and fc* is referred to k — Ha, the value when the universe scale crosses 
the Hubble horizon during inflation. 

Combining WMAP three-year data[l2] with the Sloan Digital Sky Survey (SDSS) large 



scale structure surveys [2l|, it is found an upper bound for R given by R(k* ~ 0.002 Mpc 



< 



0.28 (95%CL), where k* ~0.002 Mpc -1 corresponds to / = r k ~ 30, with the distance to the 
decoupling surface r = 14,400 Mpc. The SDSS measures galaxy distributions at red-shifts 
a ~ 0.1 and probes k in the range 0.016 h Mpc _1 < k <0.011 h Mpc" 1 . The recent WMAP 
three-year results give the values for the scalar curvature spectrum Pn{h*) = 255j I (k^)/A ~ 
2.3 x 10~ 9 and the scalar-tensor ratio R{k 1r ) = 0.095. We will make use of these values to 
set constrains on the parameters for our model. 
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V. EXPONENTIAL POTENTIAL IN THE HIGH DISSIPATION APPROACH 



Let us consider an inflaton scalar field on the brane with exponential potential. This 
potential occur naturally in some fundamental theories such as string/M theories and it is 
intensively studied in DGP models. 

We write for the exponential potential as V = VoeT^^^--, where Vo and p are two 
constants, and V{<p) — > as — ► oo, which means p > 0. An estimation of these 
parameters are give for cool inflation on a DGP brane in Ref. 10]. In the following, we will 
restrict ourselves to the high dissipation regimen i.e. r ^> 1. 



A. T oc [v((f)) + 1 + e{Al + 2r? v^)) 1 ' 2 )] 1/2 case. 



This choice is motivated by the fact that r becomes a constant. By using the exponential 
potential, we find that the slow-roll parameters become 



pr 



1 + eri(A 2 + 2rj v) 



-1/2 



[l+ 1 -(l + e(A 2 + 2r ]V y 



/2 



, and a 



pro 



l + e(^ + 2777j)V2)]^ 



pa ' 



where we have define v 

By integrating Eq. (1191) the number of e-folds results in iV = ^ [h(vf) — h(v*)], where 



h(v) = - 



In r - , - ' - ( Wv- 1 + ^-tanh.- 1 (W/A ) 

Ao 



(28) 



and W = W(v) = {Al + 2r,v) 1 / 2 . 

In section ILTTI we specified that the inflationary phase can exit naturally without any other 
mechanism in the branch e = — 1. In this case the condition for inflation end, £ — 1, gives a 
quintuple equation that can not be solved analytically. Note that the same happens in cool 
inflation for the DGP brane [lo| . 

From Eq. (12211 . when r ^> 1, we obtain that the scalar power spectrum becomes 



Pn{k) 



5/2 

Pro 
4n 2 



Po 



l /2 T r v 1 ' 2 



l + -(l + e WO 



5/2 



(29) 



and from Eq. (1271) the tensor-scalar ratio, is given by 



R(k) 



16tt 



1/2 1/2 

Po v 1 

5/2 rp 



1 -3/2 



l + -(l + eW) 

v 



coth 



2T 



(30) 



By using the WMAP three year data where Pn(k*) ~ 2.3 x 1CT 9 , R(k*) = 0.095, and 
choosing the parameter T ~ T r , we obtained from Eqs. fl29|) and (1301 that u* = 2 ^Tj 1 ^ , 
where the C is given by C ~ 5,1 x , and v u = V*/u 4 . Recall that A = =>- A) = 1 
and e = — 1. 

Now we consider the special case in which we fixe N = 60, p = 5 and ro = 10. In this 
special case we obtained that Vf = 0.05, v* = 36 and = 7.5 x 10~ 13 . Also, we choose 
the parameters K = 0.002 Mpc" 1 and T ~ T r ~ 0.24 x 10 16 GeV. We also take r? = 0.99, 
although the model can be used for any value of 77. 

From Eqs.Q and (jlj) we obtained that the tension of the brane, A becomes = ^ (1 — 

r -1 1/6 

77) ~ 2.1x10 16 , and the 5-dimensional gravitational constant results to be ^ = ~ 
0.0039. Analogously, we get that the cosmological constant 5 A in the bulk can be evaluate 
and it becomes 5 A//i 2 ~ 2.1 x 10~ 18 . We note that the A, 772,5 and 5 A parameters become 
decrease by two, one and fourth orders of magnitude, respectively, when are compared with 
those analogous cool article related to inflation in the DGP brane. 



B. T oc [v(cf>) + 1 + e(A 2 + 2r? t#)) 1/2 )] 3/2 case. 



This choice allows us to get r as a function of H in analogy with that one used by authors 



of Refs. 



22 



23j. Therefore, we take r to be of the form, r = (3 fi 2 /po)H 2 , and thus we find 



for the slow-roll parameters 



1 

pv 



1 + erj{Al + 2r) v) 



-1/2 



[1 + 1(1 + e (Al + 2rjv)y 2 
The number of e-folds is given by A" 
defined by 



and a 



2 

pv 



[v + l + e(A 2 + 2r l v) 1 / 2 )Y 



p [c> 



^s(vf) — where the function Q(v) is 



$S(v) = - 



l + A 2 e 2 



(eW 


2' 






4- - 


+ 




V 





4 ' ( ^, + " ) t M h-(H7A 



A 2 



The A" parameter has to have an appropriated values (60 or so) in order to solve the 
standard cosmological puzzles. In order to do so we need the following inequality Q^ii*) < 
$$(v f) — 120/p with Q{vf) > 120/p to be satisfied. Again, the condition for inflation becomes 
to an end is e(vf) = 1, has roots that can not be written down in closed algebraic form. 

At the epoch when the cosmological scale exit the horizon, the dissipation parameter 
becomes r(0*) = r* ^> 1 =>- p <C 3p 2 H 2 , resulting in the requirement that [v* + 1 + 
e(A 2 + 2r l v*) 1 / 2 ] > 1. 
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From Eq. (l22p . and for r > 1, we obtain that the scalar power spectrum becomes 

5 



V 
4tt 2 



pl /2 T r v* 

M 3 



l + -(l + eW) 



(31) 



k=k, 



and from Eg. (1271) the tensor-scalar ratio, is given by 



R(k) 



24 vr 
P 



1/2 

Po 
pT r v 2 



l + -(l + eW) 

v 



-4 


" k ' 




coth 






2T_ 





(32) 



By using the WMAP three year data, where P n (k*) ^ 2.3 x 10~ 9 , R(k t ) = 0.095, and 
choosing the parameter T ~ T r , we obtained from Eqs. (l3~TT) and (1321 that v* = 2 ^c-i)^ , 
where now C is given by C ~ Llxl( ! J r 1 , and ^ = K/p 4 - Here, once again we have taken 



A = =>■ A = 1 and e = -1. 

Let us consider some numerical situations. For A = 60 and p = 60, we obtain that 
Vf = 0.63, v * = 8 and f M = 2.1 x 10 -13 , where again we have taken the parameters fc* = 0.002 
Mpc" 1 , T ~ T r ~ 0.24 x 10 16 GeV and 77 = 0.99. 



From Eqs.flSD an d H), the tension of the brane becomes 



1 A 



[l-r))~ 2.6 x 10- 16 , 
~ 0.0041. 



v PO 
Gfi 4 



and the 5-dimensional gravitational constant can be written as — : 
Analogously, we obtain that 5 A//j 2 ~ 2.7 x 10~ 18 . Note that theses values are similar to 
their analogous for the case in which r = r$. 

In Fig.([T|) we plot the scalar spectrum index n s and the running spectral index a s versus 
the ratio v = V/po- In doing this, we have taken three different values of the parameter p. 
Note that in considering the n s versus v = V/po plot the WMAP-three data favors hight 
values of this parameter. However, from the plot a s versus v = V/po we could not say the 
same, since it seems that any value of p would give the same value of a s for hight enough 
v = V/po- 

Another interesting case to study is to consider T = r = constant i.e. r oc H~ l . 
Analogously as in the other cases, we take the special case with fixed N = 60 and p = 80. 
For these values, we have Vf = 0.05, t>* = 50, = 1.5 x 10~ 13 and T = 8 x 10 15 GeV. 
Here, again we have taken the parameters k* = 0.002 Mpc -1 , T ~ T r ~ 0.24 x 10 16 GeV 
and 77 = 0.99. From Eqs.(|3]) and (BJ, the tension result to be \ ~ 2.9 x 10~ 17 , and 



the 5-dimensional gravitational constant becomes — ~ 



0.0028. Finally, we obtain that 
5 A/p 2 ~ 3.0 x 10 -19 . Note that these values are similar to those values in which T is 
variable. 
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II 



1 1 1 1 1 1 1 


1 ; p 


= 80 


: ; P 


= 60 


1 ; P 

1 1 . : i . i 


= 40 



v = V /p o 

FIG. 1: The upper panel show the evolution of the scalar spectrum index n s versus v = V/po, and 
the lower panel show the evolution of the running scalar spectrum index, a s , versus v = V/pq. In 
both panels we have used different values for the parameter p, for r oc H 2 . Here, we have taken 
the values rj = 0.99, .4.0 = 1, and e = —1. 

VI. CONCLUSIONS 



In this paper we have investigated the warm inflationary scenario on a warped DGP 
brane. In the slow-roll approximation we have found a general relationship between the 
radiation and scalar field energy densities. This has led us to a general criterium for warm 
inflation in DGP brane to occur (see Eq. fll7l) ). 

Our specific models are described by an exponential potential and we have con- 
sider different cases for the dissipation coefficient, V. In the first case, we took V oc 
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1 /2 

[v((j)) + 1 + e(Al + 2r/u(0)) 1 / 2 )] . Here, we have found that the condition for infla- 
tion to end presents the same characteristic that occurs in cool inflation for the DGP 
brane [h]], except that it depends on the extra parameter Tq. For the case in which 
the dissipation coefficient T is taken to be a function of the scalar field, i.e. T oc 
[v(4>) + 1 + e(A% + 2r^f (0)) 1 / 2 )] 3//2 , it was possible to describe an appropriate warm infla- 
tionary universe model on DGP brane. In these cases, we have obtained the explicit expres- 
sions for the corresponding scalar spectrum index and the running of the scalar spectrum 
index. We have also study the situation in which r = Tq = Const. Here, we have found 
that the values are similar to those found when T is a function of the scalar field. 

By using the WMAP three year data and consider a special case with fixed N, p and 
r , we have found the values of the parameters A, 5 A and m 5 . Using the above parameters 
we will check with a numerical example that the evolution of the universe really undergoes 
a 4-dimensional stage, then a 5-dimensional stage, and finally a 4-dimensional stage again. 
From Eqs.flE]) to (J7J) the dimensional transition occurs if Ap 2 /6r??| <C 1 [ij. In our models 
that we have worked out we have obtained that A//i 4 ~ 2 x 10~ 16 , vn^j p, ~ 0.004 and 
therefore Ap 2 /6r??| = 0.008 <C 1. Also, t>* ~ p*/po ~ 10 for these models and thus the 
universe inflates in a 4-dimensional stage when the cosmic scale crossed the Hubble horizon 
during inflation, since p* ^> po ^> m \- 111 the intermediate energy region (5-dimensional), it 
is necessary that m\ <C p <C po an d using the final value of Vf ~ Pf/po, we find that pf <C 1 
and pf/m\ ~ 10. Therefore, the inflationary phase exists in a 5-dimensional stage. Finally, 
the universe stops inflating and enters in a radiation dominated Big-Bang phase, in which, 
the energy density decreases very fast and the universe becomes 4-dimensional again. 

We should note that other properties of this model deserves further study. For example, 
we did not study bispectrum of density perturbations. While cool inflation typically predicts 
a nearly vanishing bispectrum, and hence a small (just a few per cent) deviation from 
Gaussianity in density fluctuations -see e.g. 24], warm inflation clearly predicts a non- 
vanishing bispectrum. We left this question for a subsequent study of the warm-DGP brane- 
world. Also, a more accurate calculation for the density perturbation would be necessary in 
order to check the validity of expression (1201) . We intend to return to this point in the near 
future by working an approach analogous to that followed in Refs. 22J and 231 ] . 
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